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ABSTRACT 

The emission measure distribution in the upper transition region and corona of e Eri is 
derived from observed emission Hne fluxes. Theoretical emission measure distributions 
are calculated assuming that the radiation losses are balanced by the net conductive 
flux. We discuss how the area factor of the emitting regions as a function of tempera- 
ture can be derived from a comparison between these emission measure distributions. 
It is found that the filling factor varies from ^ 0.2 in the mid transition region to ~ 1.0 
in the inner corona. The sensitivity of these results to the adopted ion fractions, the 
iron abundance and other parameters is discussed. The area factors found are qual- 
itatively similar to the observed structure of the solar atmosphere, and can be used 
to constrain two-component models of the chromosphere. Given further observations, 
the method could be applied to investigate the trends in filling factors with indicators 
of stellar activity. 
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1 INTRODUCTION 

It is well known from direct observations that the solar chro- 
mosphere and transition region are far from homogeneous 
(see, for example, the uv images shown by Reeves 1976 and 
Gallagher et al. 1998). In spatially resolved observations of 
the Sun, transition region line emission tends to be dom- 
inated by contributions from the network boundaries. The 
solar-stellar analogy tells us that a similar effect must be ex- 
pected on the surfaces of other main-sequence stars. There- 
fore, it is of interest to determine the filling factor of emit- 
ting material in solar-like stars and look for correlations with 
stellar parameters, such as has been done for surface (pho- 
tospheric) magnetic filling factors by Montesinos & Jordan 
(1993). Such correlations may help to identify the physical 
processes which control coronal/chromospheric activity. 

The current generation of high resolution spectro- 
graphic instruments are providing spectra of unprecedented 
quality in the ultraviolet (uv) and X-ray wavelength ranges. 
Unfortunately, such data cannot provide direct information 
about the spatial distribution of emitting material in main- 
sequence star atmospheres since it is beyond the scope of 
current instruments to resolve dwarf star surfaces, with the 
exception of the Sun. There is little prospect of overcoming 
this limitation in the near future. 

In the absence of spatial information, semi-empirical 
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modelling of main-sequence star spectra is effectively limited 
to the creation of "mean" 1-D atmospheric models whose 
parameters are those of a hypothetical homogeneous atmo- 
sphere which reproduces the observed spectrum (e.g. the 
e Eri models constructed by Kelch 1978 or Sim 2002). Al- 
though such models are useful, tighter constraints on the 
true atmospheric parameters, and therefore the nature of the 
physical processes which operate, could be achieved if the 
distribution of emitting material were known. (Information 
on the structure within the atmospheres of cool giant stars 
has been obtained through studies of molecular fiuorescence, 
see McMurry & Jordan 2000.) The approach which could 
be taken is illustrated by the work of Cuntz et al. (1999). 
They have explored two-component chromospheres of K2 V 
stars, in which one component is heated by longitudinal tube 
waves and the other by acoustic waves. They simulate the 
supergranulation network using a fiux tube model based on 
empirical correlations between observed surface magnetic 
fiuxes and stellar rotation rates (Cuntz, Ulmschneider & 
Musielak, 1998) and calculate the flux in the Ca ll H and 
K lines expected from magnetoacoustic and acoustic waves. 
In the absence of observational constraints they use an av- 
erage chromospheric network filling factor of 0.305. Thus a 
method which gives the filling factor in the transition region 
between the chromosphere and corona would provide valu- 
able constraints on the above type of theoretical modelling. 

Kopp (1972) discussed the energy balance in the solar 
transition region, in the light of early observations of the 
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solar supergranulation structure. He introduced an area fill- 
ing factor and by comparing emission measure distributions 
based on spatially integrated observations with the predic- 
tions of the energy balance equation he proposed that the 
transition region emission is restricted to an area of about 
20 per cent of the solar surface area. In this case, the con- 
ductive flux in the mid-transition region is reduced and the 
lower transition region requires an additional source of heat- 
ing. This type of structure was also explored by Gabriel 
(1976). The influence of the cross-section area and magnetic 
field expansion factor has also been recognized in the con- 
text of loop models of solar and stellar active regions (see 
Underwood, Antiochos & Vesecky 1981). 

A method of deriving area factors has also been dis- 
cussed by Jordan (2000), in which the total pressure is not 
kept constant but is allowed to vary according to hydrostatic 
equilibrium. Although Philippides (1996) used this approach 
to estimate the filling factor at Te — 10^ K in several main- 
sequence stars, the observational data available were not 
sufficient to find the variation of the flUing factor with Te- 
Here we apply the method to spatially unresolved stellar ob- 
servations to derive, for the first time, the area filling factor 
as a function of Te in the upper transition region of e Eri. 

The star selected for this analysis, e Eri (K2V), is a 
nearby dwarf that has moderately high levels of chromo- 
spheric/coronal activity and for which high quality spectro- 
scopic data are available. This paper forms part of our study 
of the outer atmosphere of e Eri: see Jordan et al. (20Qla), 
Jordan et al. (2QQlb), Sim & Jordan (2003) and Sim (2003, 
in preparation) for details. 

In the next section the theory required for the analysis 
is discussed. Section 3 discusses the observations that are 
used to constrain the emission measure distribution (Sec- 
tion 4). Section 5 presents "first-cut" area filling factors de- 
duced from the data and discusses possible sources of error. 
In Section 6, a self-consistent determination of the filling 
factors is presented. The results are discussed in Section 7 
and conclusions are drawn in Section 8. 



2 THEORY 

Spatially unresolved stellar observations can be used to con- 
strain the apparent mean emission measure (see Pan & Jor- 
dan 1995) distribution defined by 



Em°-'{n) = /^"'^ N4T)NH{T)^G{r)f{r)^ dT 

(1) 

where Te is the electron temperature at radius r, and 
Nh are the electron and hydrogen number densities respec- 
tively, /(r) = / Ri (i?, is the photospheric stellar radius 
and /(r) = 1 at the stellar surface). The dilution factor 
G(r) specifies the fraction of photons (emitted in an opti- 
cally thin line at radial coordinate r) that escape the star 
uninterrupted by the stellar surface. For spherical geometry. 



G(r)= 0.5 (l-yi-[,fW]-') 



(2) 



The 0.3 included in the notation for the emission measure in 
equation (1) is to remind the reader that the emission mea- 



sure adopted is an integral over a logarithmic temperature 
range AlogTe = 0.3 dex. This temperature range is chosen 
since it is the typical temperature range across which in- 
dividual spectral lines form in the transition regions of cool 
stars. The area filling factor A(r) /A, (r) specifies the fraction 
of the total stellar surface area at radius r (A*(r) = A-kt'^) 
which is occupied by emitting material. The purpose of this 
paper is to show how values for A[r) /A^ (r) can be extracted 
from spatially unresolved data. 

In the transition region where the pressure varies slowly 
as a function of the temperature, equation (1) can be written 
approximately as 



G{r)f(r] 



A(r) dr PePff 
A,{r) dTe ^klTe 



(3) 



where Pe = NeknTe and Ph ~ NnksTe are the electron 
and hydrogen pressures respectively, and are approximately 
constant over the region where an individual line is formed. 
The approximation that dr/dTe is constant in the region 
of line formation is less good, but is improved upon when 
finding final models by iteration. 

As discussed by Jordan (2000), it is possible to explain 
the complete radiation losses from the upper transition re- 
gion in terms of the divergence of the net conductive energy 
flux from the overlying corona. The justification for this is 
that the emission lines are formed over a pressure-squared 
scale height ~ 40, 000 km, whereas the main heating of the 
quiet corona occurs at much greater heights. It cannot be 
ruled out that there is an additional source of heating in this 
part of the atmosphere, but it is certain that any additional 
heating must be small. The conductive flux can be expressed 
as 



5/2 



dTe 

dr 



(4) 



where «: « 1.1 x 10"® ergs cm'^ s"^ K'^/^ _ ;^o^ K 

(Spitzer 1956). k increases at lower temperature (by about 
~ 30 per cent by 2 x 10^ K), but it is taken as constant here 
for simplicity. Building on the simpler plane parallel approx- 
imation used by Jordan & Brown (1981), the logarithm of 
equation (3) can be differentiated and the temperature gra- 
dient replaced by the conductive ffux using equation (4) . By 
assuming an energy balance between conduction and radia- 
tion, the net conductive ffux can be replaced by the radiation 
losses and the following equation obtained: 



dlogT, 



log 



A,(r)£mO 



+ 



+ 2 



Mr) 
dlogyl(r)/^,(r) 
dlogTe 

2Prad(re) 

KP,PHTi'''G{rYf{rY 



dIogG'(r 
dlogTe 



+ 2 



dlog/(r) 
dlogTe 



dlogPg 
dlogTe 

A,{r)Em° 



+ 2 



Mr) 



(5) 



Trad(Te) IS the radiative power loss function which describes 
the emissivity of a plasma of fixed composition at tempera- 
ture Te. Pg = NgksTg is the total gas pressure. 

This equation can be solved numerically to obtain the 
distribution of A,{r)EmP'''^ / Air) as a function of Te. By 
comparison with values of EniP'''^{Te) (deduced from obser- 
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vations) this allows the area filling factor A{r) /Ai,{r) to be 
derived. 



3 OBSERVATIONS 

The mean emission measure distribution can be constrained 
using observed emission line fluxes. The details of this pro- 
cedure are well known and have been discussed previously 
by various authors (e.g. Jordan & Brown 1981, Jordan et al. 
1987; Griffiths & Jordan 1998). 

In this paper, the mean emission measure distribution in 
the upper transition region of e Eri is constrained using the 
observed fluxes of emission lines of iron. For the most part, 
the observed fluxes used are from the analysis performed 
by Schmitt et al. (1996) using data recorded by the Extreme 
Ultraviolet Explorer (EUVE) satellite. The data and line flux 
measurements are discussed by Schmitt et al. (1996) and 
Laming, Drake & Widing (1996). Schmitt et al. (1996) give 
fluxes for Fe lines from ionisation stages ix - xvi and xvill - 
XXI. For the current work, only one line from each ionisation 
stage detected in the EUVE data is used. For Fe xil the 
three lines at around 193 A have been added to give the 
total flux in the *S -*P multiplet. In other cases where there 
are several lines from the same ion reported by Schmitt et 
al. (1996), the strongest line has generally been used. The 
lines used are listed in Table 1 together with the percentage 
errors in the measured fluxes given by Schmitt et al. (1996). 

The observed EUVE fluxes have been converted to 
fluxes at the stellar surface using the stellar parameters for 
e Eri listed in table 1 of Jordan et al. (2001a). Subsequent to 
the work carried out by Schmitt et al. (1996), a more accu- 
rate determination of the hydrogen column density along the 
line of sight to e Eri has been made by analysis of interstellar 
absorption features in the Lyman a line profile (Dring et al. 
1997). This gives a column density of log iV(cm"^) = 17.875, 
which is smaller than that used by Schmitt et al. (1996) by 
close to a factor of 1.5. The wavelength dependent transmis- 
sion of the interstellar medium with the new column density 
has been taken from calculations by Philippides (1996) and 
has been used to correct the observed fluxes for interstellar 
absorption. 

The EUVE data are supplemented by our recent detec- 
tion of the forbidden ultraviolet Fe xil lines in STIS spectra 
of e Eri. These data, and the measurement of the Fe xil line 
fluxes, have been presented elsewhere (Jordan et al. 2001a). 



4 EMISSION MEASURE LOCI 

As discussed by Jordan & Brown (1981), each line flux can 
be used to construct a locus which places an approximate 
upper bound on the mean emission measure distribution: 



(6) 



where is the stellar surface flux for the line under con- 
sideration and K(Te) is a contribution function defined by 
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Figure 1. Apparent emission measure loci deduced using the ion- 
isation balance of Arnaud & Rothenfiug (1985). The solid curves 
were computed using the observed EUVE fluxes of Schmitt et al. 
(1996) and are labelled with the ionisation stage of iron which 
they represent. The dashed curve was computed using the STIS 
Fe xil 1242-A flux reported by Jordan et al. (2001a). 



where A is the wavelength of the line, Aui is the Einstein 
coefficient for spontaneous emission in the line, and the A'''s 
are number densities with the subscripts u, I and E refer- 
ring to ions in the upper level in the transition, the ionisa- 
tion stage in which the transition occurs and the element 
giving rise to the transition respectively. K{Te) has been 
calculated for each of the lines listed in Table 1, and for 
the Fe xil 1242-A line. The CHIANTI atomic database (v4; 
Dere et al. 1997; Young et al. 2003) was used to compute 
Nu/Ni as a function of Te.The calculations were performed 
using the best value of the mean transition region electron 
pressure logPe(cm~'' K) = 15.68 obtained by Jordan et al. 
(2001b). Two sets of calculations for the ionisation fractions 
of iron have been used (Arnaud & Rothenfiug 1985; Arnaud 
& Raymond 1992) and the stellar photospheric iron abun- 
dance [Fe/H]=-0.09 found by Drake & Smith (1993) has 
been adopted. (The solar iron abundance was taken from 
Grevesse & Sauval 1998). Fig. 1 shows the loci F^/K{Te) 
derived from the observed fiuxes using the ionisation bal- 
ance calculations of Arnaud & Rothenfiug (1985) and Fig. 2 
shows the loci obtained with the ionisation balance of Ar- 
naud & Raymond (1992). 

As expected, the emission measure loci suggest a 
smooth emission measure distribution which increases with 
temperature up to a peak at several million degrees. The 
temperature of the peak emission measure is difficult to 
determine since there are no Fe xvil lines in the data 
which would span the critical temperatures range. Compar- 
ing Figs. 1 and 2 shows that the choice of ionisation balance 
affects the form of the emission measure distribution. In par- 
ticular there are significant differences in the ion populations 
of Fe IX, XVI and xvill which alter the apparent shape. The 
Arnaud & Raymond (1992) ionisation balance calculations 
suggest a higher peak temperature (logTe(K) ~ 6.6) than 
those of Arnaud & Rothenfiug (1985) (logTe(K) ~ 6.4). 
In their analysis. Laming et al. (1996) adopted the Arnaud 
& Rothenfiug (1985) ionisation balance (leading them to a 
coronal temperature logTe(K) ~ 6.4) since solar data sug- 
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Table 1. Lines of iron observed in EUVE spectra of e Eri. 



Ion 


Ao (A)- 


Transition 


% Error* 


Fe IX 


171.073 


^So - 3p^3d 


17 


Fe X 


174.534 


3s23p5 2p3/2 - 3s23p*(3p)3d ^Dg/a 


18 


Fe XI 


180.408 


3s23p^ ■^Pa - 3s23p-^(^S)3d •^Ds 


19 


Fe XII 


192.393/193.521/195.118 


3s23p3 433/2 - 3s23p2(3p)3d ''Pi/2,3/2,5/2 


40 


Fe XIII 


203.828 


3s23p^ ■^Pa - 3s23p3d ^Ds 


38 


Fe XIV 


211.317 


3s-^3d ^Pi /o - 3s^3d ^Do/o 


18 


Fe XV 


284.160 


3s2 iSo - 3s3p iPi 


7 


Fe XVI 


335.410 


3s 2Si/2 - 3p 2P3/2 


11 


Fe XVIII 


93.923 


2s22p5 2P3/2 - 2s2p6 2Si/2 


22 


Fe XIX 


108.356 


2s22p4 3p2 - 2s2p5 3p2 


24 


Fe XX 


132.841 


2s22p3 433/2 - 2s2p4 4pg/2 


21 


Fe XXI 


128.752 


2s22p2 3po _ 2s2p3 ^Di 


40 


''The error 


column gives the percenta 


" Rest wavelength. 

ge error in the measured fiux from Schmitt et al. (1996). 
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Figure 2. The same as for Fig. 1, but using the ionisation balance 
of Arnaud & Raymond (1992). 



EUVE Fe XII loci shown in Figs. 1 and 2 are computed using 
the fluxes for all three lines in the *P - multiplet from 
Schmitt et al. (1996). If the three fluxes reported by Schmitt 
et al. (1996) are used separately, the 195-A line gives a locus 
that is significantly higher than that obtained by combining 
all three lines of the multiplet while the 192-A line gives 
a locus that is significantly lower. The 193-A line gives a 
locus that agrees well with that of the complete multiplet. 
It is plausible that the discrepancy (factor of 1.5) between 
the EUVE and STIS Fe xil loci is the result of systematic 
errors in the atomic data, but the measurement errors for 
the EUVE fluxes (see Table 1) are large enough to account 
for most of the difference. As suggested by Jordan et al. 
(2001a), there may be a real difference related to variations 
in coronal activity between the times at which the two sets 
of observations were made, but the uncertainty in the EUVE 
fiux is too large to establish this with certainty. 



gests that these calculations are more appropriate than those 
of Arnaud & Raymond (1992) for Fe ix (see the discussion 
by Laming, Drake & Widing 1995) . In this paper, all calcu- 
lations are carried out with both sets of ionisation fractions, 
and the differences between the results are assumed to be 
indicative of the errors in the ionisation balance. 

Above the peak temperature the emission measure 
drops with increasing temperature. There appear to be some 
discrepancies between the loci for the high temperature ions. 
In particular the Fe xx locus appears to lie too high com- 
pared to those of Fe xix and Fe xxi. Material at tempera- 
tures which are significantly greater than that of the peak 
emission measure most likely occurs in stellar active regions, 
where it is magnetically confined in hot loops (such as dis- 
cussed by Schmitt et al. 1996) and will not be discussed 
further here. Solar observations support this interpretation. 

The locus from the 1242-A STIS line hes above that 
from the EUVE Fe xil lines by about a factor of 1.5. This is 
significant but much smaller than the discrepancy reported 
by Jordan et al. (2001a) (they reported a factor of 3 dif- 
ference between the STIS and iJ 17 results). Jordan et al. 
(2001a) computed the EUVE Fe xii locus using the fiux re- 
ported by Laming et al. (1996) for the 195-A line only. The 



5 EMISSION MEASURE MODELS 
5.1 Models 

As discussed in Section 2, equation (5) can be used to con- 
struct theoretical emission measure distributions under the 
assumption of an energy balance between thermal conduc- 
tion and radiation in the upper transition region. To de- 
termine a unique emission measure distribution from equa- 
tion (5) two boundary conditions are required: the pressure 
and emission measure must both be fixed at some tempera- 
ture. The pressure at all other points is then determined by 
enforcing hydrostatic equilibrium. For this, a contribution to 
the pressure from turbulent motions is included, calculated 
using a most-probable turbulent velocity ^ « 21 km s~^ 
which is appropriate for the temperature range under con- 
sideration (Sim & Jordan 2003). (Here, this turbulent pres- 
sure contribution is very small.) With the pressure so deter- 
mined, equation (5) can be integrated numerically to deter- 
mine the emission measure at all other points. 

The first boundary condition adopted is that the elec- 
tron pressure at logTe(K) = 5.3 should be consistent with 
the mean value determined by Jordan et al. (2001b). The ef- 
fects of slightly varying the choice of fixed electron pressure 
at logre(K) = 5.3 have been explored (see below). 
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It is a natural consequence of the signs of the terms 
in equation (5) that the computed emission measure distri- 
bution will pass through a minimum in the transition re- 
gion when the magnitude of the last term in equation (5) 
is close to 3/2. Therefore, the second boundary condition 
has been chosen to be that this minimum should occur at 
logre(K) = 5.3, which is approximately the temperature at 
which the minimum in the emission measure distribution is 
observed to occur in the Sun (see e.g. Macpherson & Jor- 
dan 1999). The emission measure distribution for e Eri (Sim 
2002) shows a decrease up to at least logr5(K) = 5.3. 

Using these two boundary conditions, emission measure 
distributions {At{r)Em'^'^ / A{r)) have been computed using 
equation (5). At this stage, the third term on the right- 
hand side (RHS) of equation (5) which accounts for the 
variation in the area factor with Te has been ignored. The 
effect of this term is discussed in Section 6. The distribu- 
tions have been calculated from an upper boundary temper- 
ature logTmax ~ 6.5 down to a lower boundary temperature 
logTmin = 5.3. Distributions of A,{r)Em°-^ /A{r) were cal- 
culated using a range of starting values for A* {r)Em'''^ /A(r) 
and the electron pressure at the top point in the model 
in order to find the distributions which satisfies the chosen 
boundary conditions (see above). Two of the computed dis- 
tributions are shown in Fig. 3 (the solid and dotted lines). 
These distributions differ only in the choice of fixed pres- 
sure at the bottom end. Fig. 4 shows the electron pres- 
sure (Pe = NeTe) as a function of Te for the same mod- 
els. The first model (solid line in the figures) has a pressure 
which is slightly lower than, but close to, the best pressure 
(logPe(cm~'^ K) ~ 15.68) while the other (dotted line) has 
a higher pressure, just within the error margins of the best 
pressure (logPe(cm-^ K) ~ 15.78). 

In both models, it can be seen that above the mid- 
transition region the theoretical emission measure increases 
roughly proportional to T^^^^ owing to the dominance of the 
penultimate term in equation (5). At higher temperatures, 
the emission measure begins to flatten off and eventually 
starts to decrease. This is the result of the negative contri- 
butions to the RHS of equation (5) from the 1st and 4th 
terms which become significant at high temperatures (the 
last term is small at high temperatures). 

5.2 Filling factors 

The theoretical distributions plotted in Fig. 3 can be used 
to estimate fractional filling factors for each of the observed 
iron lines. The apparent emission measure distribution can 
be used to compute approximate line fluxes via 

F*^ J K{n)Em'^ ''(n) d\ogT,/0.3 (8) 

where the integral runs over the whole temperature range 
of the model. If the apparent emission measure in equa- 
tion (8) is replaced with the computed distribution of 
At,{r)Em°-^ /A{r) then P«4*(r)/4(r) is obtained instead of 
P*. This can be compared to the observed value of the line 
flux to deduce A(r)/A*(r) in the region of line formation. 
The value so deduced is only an estimate because of the ap- 
proximate nature of equation (8) and must be interpreted as 
a weighted average of the area filling factor over the region 
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Figure 3. Theoretical curves of Af(r")Em^''^ j Air) against tem- 
perature. The dash-triple-dot curve is the Fe XII STIS locus from 
Fig. 1. The solid and dotted curves are for models in which the 
variation of the fiUing factor with temperature (on the RHS of 
eqn. (5)) has not been included, and refer to different (higher 
and lower) pressures at logTe = 5.3. The dashed and dashed- 
dot curves are for the corresponding self-consistent models which 
include the variation of the filling factor with temperature (see 
Section 6). 



1 6.0 ^ 

15,8 - 




log Te (K) 

Figure 4. The electron pressure as a function of the electron 
temperature for the models plotted in Fig. 3, using the same 
symbols. 

of line formation. It should be a good estimate provided 
that the physical extent of the region of line formation is 
reasonably small (as is the case for all the lines considered 
here) . 

The filling factors deduced from the distributions dis- 
cussed above are given in Table 2. Filling factors are only 
given for individual lines that form in the upper transition 
region (Fe ix - Fe xvi) since the higher temperature lines 
probably form under different conditions in hot active re- 
gion loops (see Section 4). For each line the approximate 
temperature of line formation (determined from the emis- 
sion measure loci) and four filling factors are given. These 
correspond to the four possible combinations of the two 
models and the two sets of ionisation balance calculations 
(Arnaud & Rothenfiug 1985; Arnaud & Raymond 1992). In 
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Figure 5. Area filling factors deduced from the model with 
logPe ~ 15.78 using the Arnaud & Raymond (1992) ionisation 
balance calculations. The open diamonds show the results from 
the EUVE iron lines. The open triangles show the results from 
the STIS data (the Fe XII 1242-A line and the mid-transition re- 
gion lines). These results do not include the term in equation (5) 
describing the variation of the area factor with temperature. The 
filled symbols show the results when the variation of the filling 
factor with temperature is included in a self-consistent manner 
(see Section 6). 

addition to the filling factors deduced from the upper tran- 
sition region lines, filling factors are given in Table 2 for the 
mid-transition region (logre(K) ~ 5.3) which were deduced 
by comparing the computed Ai.{r)EmP'^ / A{r) with values 
of EmP'^ from a mean transition region emission measure 
distribution constructed from STIS data (Sim 2002). This 
distribution will be discussed in a forthcoming paper (Sim 
2003, in preparation). 

The models suggest that the mid transition region fill- 
ing factor is several tens per cent and that the filling factor 
starts increasing with temperature above logre(K) ~ 5.9, 
becoming close to 1.0 at coronal temperatures. To illustrate 
this point. Fig. 5 shows (open symbols) a plot of the area 
filling factors deduced from the model with higher pressure 
using the Arnaud & Raymond (1992) ionisation balance (the 
error bars indicate the random measurement errors but do 
not attempt to address the potential systematic errors which 
are discussed below). This trend with temperature is strik- 
ing and consistent with what would be expected from the 
standard model of expanding magnetic funnels in the upper 
atmosphere (Gabriel 1976). 

The filling factors are not strongly sensitive to the 
adopted boundary pressure (they are 20 - 40 per cent 
smaller at the larger pressure). As discussed above, how- 
ever, several of the filling factors are quite sensitive to the 
adopted ionisation balance. 

On physical grounds, the filling factor cannot be greater 
than 1.0. In all cases the filling factor calculated for Fe xvi 
is greater than 1.0 and, for all cases except the high pres- 
sure/Arnaud & Raymond (1992) combination, the filling fac- 
tor is unphysical for some of the other lines. In most cases 
the random measurement errors are not sufficient to explain 
these unphysical results, and they point towards a system- 
atic error in the calculations. Several possible sources of sys- 
tematic errors are discussed below. 



5.3 Systematic errors 

5.3.1 Ionisation balance 

The quality of the atomic data is probably the most signif- 
icant source of uncertainty, particularly those used in the 
ionisation balance. By comparing the filling factors deduced 
using Arnaud & Rothenflug (1985) and Arnaud & Raymond 
(1992) it is clear that for several of the ions (Fe ix, Fe xv 
and Fe xvi) the uncertainty in the ion balance may be com- 
parable to, or even greater than, a factor of 2. For Fe xil 
and Fe xill the uncertainty is smaller (~ 40 per cent) and it 
is negligible for Fe x and Fe xi. 

5.3.2 Hydrogen column density 

The results are sensitive to the adopted interstellar hydrogen 
column density since it determines the wavelength depen- 
dent transmission factor of the interstellar medium. This is 
most important for Fe xvi (the longest wavelength line). The 
error on the column density quoted by Dring et al. (1997) 
is only 20 per cent, however, and so this should not be a 
dominant contribution to the uncertainty. 

5.3.3 The abundance of iron 

A potentially important source of systematic error is the 
adopted iron abundance. The quoted error in the photo- 
spheric abundance derived by Drake & Smith (1993) is only 
0.05 dex, but it is possible that the iron abundance in the 
upper transition region is different from that in the photo- 
sphere. It has been suggested that the abundances of ele- 
ments with low (< 10 eV) first ionisation potentials (FIPs) 
are enhanced in the outer parts of cool star atmospheres 
relative to elements with high FIPs (Crawford, Price & Sul- 
livan 1972 first identified this effect in the solar wind). This 
is known as the FlP-effect (see Jordan et al. 1998 and refer- 
ences therein for a full discussion) . Laming, Drake & Widing 
(1996) were unable to conclusively detect an enhancement 
of the iron abundance in the corona of e Eri (as would be 
expected if the FlP-effect were acting), however they could 
not rule out an enhancement of similar magnitude to that in 
the Sun (increase by a factor of 2 - 3). If the upper transition 
region/coronal abundance is enhanced relative to the photo- 
spheric value adopted in the models, then the predicted area 
filling factors in Table 2 will all be systematically overesti- 
mated. (It can be shown from equation (5) that an increase 
in the adopted iron abundance by a factor x would lead to 
a reduction in the derived area factors by approximately a 
factor of ^/x.) 

5.3.4 The stellar surface gravity 

Errors in the assumed value for the stellar surface gravity g 
will affect the calculated area filling factors in a systematic 
way. In particular, if the adopted gravity is too large, the 
area factors found for the high temperature lines will be too 
large. This is because the magnitude of the 4th term on the 
RHS of equation (5) is proportional to the adopted gravity. 
Since that term is negative, if the adopted gravity is too 
large it will mean that the peak in the calculated emission 
measure distribution will occur at too low a temperature 
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Table 2. Filling factors deduced from the observed e Eri lines and the first two models shown in Fig. 3. Columns 4-7 give the filling 
factors deduced from the different combinations of models (higher/lower transition region pressure) and ionisation balance calculations 
(Arnaud & Rothenflug 1985 (AR85) and Arnaud & Raymond 1992 (AR92)). The first part of the table gives the results from the lines 
observed with EUVE. The second part shows the results from the Fe XII 1242-A line (from STIS) and also from comparison with the 
mean emission measure distribution in the mid transition region (see text). Here logg = 4.75 is used (from Drake & Smith 1993). 
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Figure 6. The same as for the solid and dotted lines in Fig. 3, 
but using a lower gravity of logg = 4.65. 



(recall that the 4th term is important in determining where 
the peak occurs) . To investigate this effect models have been 
computed using a lower surface gravity, logg = 4.65, which 
is consistent with the error bar on the value of the surface 
gravity found by Drake & Smith (1993) (logg = 4.75 ±0.1). 
The emission measure distributions obtained with the lower 
gravity (using the same lower boundary pressures as before) 
are plotted in Fig. 6. These emission measure distributions 
have been used to compute area filling factors (using the 
Arnaud & Raymond 1992 ionisation balance) for comparison 
with those found in Table 2. The area factors deduced for 
the EUVE lines are given in Table 3. 

Comparing the area factors in Table 3 with the ap- 
propriate entries in Table 2 shows that, as expected, the 
high temperature lines have the greatest sensitivity to the 
adopted gravity: for Fe ix - Fe xil the difference in the 
adopted gravity decreases the computed area factor by less 
than 10 per cent, but for Fe xvi the reduction is close to 50 
per cent. 

Therefore, it is possible that the unphysical area factors 



for the high temperature Fe xvi line could arise because the 
adopted surface gravity is too large. (Note that earlier work, 
e.g. by Kelch 1978 did adopt the significantly lower surface 
gravity of logg = 4.5.) 



5.3.5 The radiative power loss function, PradiT^) 

In all the calculations discussed here, the fit to the radiation 
loss function 



Prad(Te) = 1.25 X IQ-^'^T, 



-1 3 - 

ergs cm s 



(9) 



derived by Philippides (1996) has been adopted. Across the 
temperature range under consideration, this fit is accurate to 
within a factor of ~ 2. For future work, in which it is hoped 
that the other sources of systematic error may be eliminated, 
it would be valuable to obtain a more accurate power loss 
function using modern atomic data. For the present, how- 
ever, this simple form is sufficiently accurate given the un- 
certainties in the ionisation fractions and iron abundance, 
and that the last term in equation (5) is only important at 
the low temperature end of the models. 



5. 3. 6 Discussion of systematic errors 

To conclude this section, it seems probable that the uncer- 
tainties in the deduced filling factors are dominated by those 
in the adopted ionisation balance calculations. The uncer- 
tainty in the Fe xvi filling factor is likely to be the largest (a 
multiplicative factor of 2 - 3) since it is the most sensitive to 
errors in the adopted hydrogen column density, surface grav- 
ity and has amongst the most uncertain ionisation fraction. 
These sources of error are certainly large enough to explain 
the moderately unphysical filling factors deduced for this 
line. 

The errors are likely to be smallest in the Fe x, Fe xi, 
Fe XII and Fe xill filling factors since these are the least 
sensitive to the systematic errors discussed above. 

The scale of the uncertainties in all cases is small enough 
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Table 3. Filling factors deduced from the iron lines observed in e Eri and the models shown in Fig. 6. Columns 4 and 5 give the filling 
factors deduced from the two models (higher/lower transition region pressure) using the ionisation balance calculations of Arnaud & 
Raymond (1992) for the EUVE lines. Here logg = 4.65 is used. 
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that the filling factors presented in Table 2 can be regarded 
as reasonable estimates. That the higher pressure model 
with the Arnaud & Raymond (1992) ionisation balance gives 
the "most physical" set of filling factors could be interpreted 
as meaning that this is the best combination of pressure and 
ionisation balance, but the errors are large enough that this 
cannot be clearly established. 

The detected trend with temperature is not expected 
to be the result of any of the errors discussed here and is 
ahnost certainly real: A{r)/At:{r) varies systematically by 
more than a factor of 5 over the temperature range consid- 
ered. This variation is larger than that expected for any of 
the uncertainties discussed above. 



6 SELF-CONSISTENT CALCULATION 

The estimates for A{r)/At{T) presented above are not self- 
consistent. This is because when equation (5) was used the 
3rd term on the RHS was neglected, but the implied filling 
factors suggest a significant variation of A(r)/A*(r) with 
logTe. To investigate this inconsistency, two more models 
have been constructed accounting for the previously ne- 
glected term. 

To do this, the area filling factor must be specified 
as a function of temperature in the modelling. In princi- 
ple self-consistency between the variation of the area factor 
adopted in the modelling and the derived area factors could 
be achieved iteratively, but since there are significant uncer- 
tainties (see above) a rigorous iterative solution is not war- 
ranted and only approximate self-consistency is pursued. To 
this end a simple form for A{r) / Af{r) as a function of Te 
has been assumed in the calculation of the two new models, 
based on the approximate area factors presented in Table 2. 
Given the form of the relevant term in equation (5), it is 
computationally convenient to express the relationship as a 
power law. The relationship adopted in the calculations is 



lofi 



Mr) 



-0.7 Iogre<5.8 

-6.5 + log 5.8 < log Te < 6.5 



(10) 



This relationship crudely describes the form suggested 
by the results in the last column of Table 2: constant below 
logTe = 5.8 and then increasing to approximately 1.0 in the 
corona. Distributions of At{r)EmP'^ / A{r) calculated from 
the two new models (which, like the first two models difi'er 



only in the adopted lower boundary pressure) are plotted 
in Fig. 3 (the dashed and dashed-dot lines). The electron 
pressures in these models are plotted in Fig. 4. 

Comparing the distributions in Fig. 3 shows the effect of 
the area variation term: it steepens the predicted emission 
measure gradient before the peak (the kink in the dashed 
and dashed-dot lines is not physical but due to the discon- 
tinuity in equation 10). Also, above the peak the emission 
measure declines more rapidly. This is because the larger 
emission measures around the peak lead to a more rapid 
decline in the pressure with temperature, making the mag- 
nitude of the 4th term on the RHS of equation (5) greater. 

Table 4 gives the area factors deduced from the new 
models in the same format as those given in Table 2. The 
area factors deduced from the new model with the higher 
pressure and the Arnaud & Raymond (1992) ionisation bal- 
ance are shown in Fig. 5 (filled symbols). A curve showing 
the assumed A(r)/Ai,{r) (equation 10) is over-plotted for 
comparison. 

Given that there are several sources of systematic error, 
the agreement between the form of A(r-) / At{r) used in the 
modelling and the values of A{r)/Af{r) predicted by the 
model is sufficiently good that the higher pressure model 
can be regarded as self-consistent if the Arnaud & Ray- 
mond (1992) ionisation balance is adopted. It is noteworthy 
that the largest discrepancies are for Fe ix, Fe xv and, in 
particular, Fe xvi which are the ions that show the largest 
uncertainties in the ionisation balance. The area factors de- 
duced from the model with lower pressure are slightly less 
self-consistent than those obtained from the higher pressure 
model, being systematically larger than those assumed in 
the modelling. Given the scale of the other sources of error, 
however, the agreement is adequate. 

Comparing Tables 2 and 4, it is clear that the main ef- 
fect of including the term describing the variation of the area 
factor with temperature is to reduce the implied area factors 
for most of the high temperature lines, leading to fewer un- 
physical {A{r) / At{r) > 1.0) values. This is the result of the 
steeper increase of emission measure with temperature. The 
only exception is Fe xvi for which the implied area factor 
is larger when the variation of the area factor is included in 
the model if the Arnaud & Raymond (1992) ionisation frac- 
tions are used. This is because Fe xvi forms at temperatures 
above the peak of the emission measure distribution, if the 
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Table 4. As for Table 2, but using the self-consistent models shown in Fig. 3. 
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Arnaud & Raymond (1992) ionisation balance calculations 
are adopted. 

In summary, the self-consistent treatment of the varia- 
tion of the area factor presented in this section indicates that 
neglecting this variation in the modelling may cause the area 
factors in the upper transition region to be systematically 
overestimated to a small, but significant, extent. Accounting 
for the variation of the area factors does not, however, alter 
the apparent trend that the area factor increases across the 
upper transition region. For now it is concluded that the 
area factors presented in Table 4 (or equivalently described 
by equation 10) are the best current estimates of the area 
filling factors in the upper transition region of e Eri and that 
most of them are accurate to within a factor of two. 



7 DISCUSSION AND PROSPECTS FOR 
FUTURE DEVELOPMENT 

In Sections 5 and 6 it has been shown that the theory de- 
scribed in Section 2 can be used to provide estimates of the 
area filling factors of emitting material at various tempera- 
tures in stellar transition regions. 

The results show (Fig. 5) that in e Eri the area filling 
factor increases with temperature markedly in the upper 
transition region, from ~ 0.2 at logTe < 5.8 to ~ 1.0 in the 
corona. There is no evidence for significant variations across 
the mid-transition region 5.3 < logTe < 5.8. 

The variation of the area factor derived for e Eri is sim- 
ilar to that found from observations of the solar supergranu- 
lation network made with the Solar and Heliospheric Obser- 
vatory {SOHO) (Gallagher et al. 1998). They find (see their 
figs. 5 and 8, in particular) that at the resolution of the Coro- 
nal Diagnostic Spectrometer (CDS) (~ 2 arcsec by 2 arcsec) 
the area occupied by the emitting material increases dra- 
matically in the upper transition, but does not vary signifi- 
cantly across the mid-transition region. These results are in 
general agreement with earlier observations and the model 
by Gabriel (1976). When the method used here is applied 
to CDS observations of the average network, a filling factor 
of ~ 0.2 at Te ~ 2 X 10^ K is deduced within the network 
(Smith & Jordan 2002), which typically occupies ~ 50 per 
cent of the solar area at this temperature (Gallagher et al. 
1998). The filling factors which we have derived are based 



on the assumption that the observed emission is dominated 
by the regions occupying the small areas. These could refer 
to structures within the supergranulation network of e Eri 
and may not be comparable to the total solar network area 
observed with the resolution of CDS. In this case, the area 
factors are larger than those found for the Sun. 

Using the solar analogy, the area factor found at Te — 
2 X lO'^ K can be applied to the lower transition region. 
This increases the intrinsic radiation losses and since ther- 
mal conduction now plays no role below ~ 10^ K, an ad- 
ditional source of heating is required. As found by Sim & 
Jordan (2003), if the observed non-thermal emission line 
widths are associated with an Alfven wave (or fast-mode 
wave) fiux, there is ample energy to heat both the corona 
and the lower transition region, although uncertainties con- 
cerning wave propagation and dissipation remain to be re- 
solved. 

The filling factor of ~ 0.2 found for the mid-transition 
region of epsilon Eri can be compared with the value of 
0.3 assumed for the chromosphere by Cuntz et al. (1999). 
Since the filling factor of magnetic flux tubes should increase 
with height, our results suggest that a smaller chromospheric 
filling factor of < 0.2 would be more appropriate; this could 
be explored in future two-component models. 

The accuracy of the results presented in this paper is 
limited by the various systematic sources of error which were 
discussed in Section 5.3. In particular the analysis of larger 
sets of lines, ideally covering a wide range of elements and 
ionisation species would help. This would reduce the sys- 
tematic errors associated with uncertainties in the ionisation 
fractions for iron and possibly help to identify possible er- 
rors arising from the adopted elemental abundance. Suitable 
data for such analysis are gradually becoming available for 
a range of stars thanks to the current generation of X-ray 
satellite missions. X-ray data are also advantageous since 
their analysis is less sensitive to the assumed interstellar 
hydrogen column density. Improved calculations of atomic 
data for iron, and the radiation loss function P^d would also 
make the modelling more quantitatively reliable. 

In future work this technique should be applied to a 
range of late-type stars with differing levels of activity. By 
doing so it should be possible to look for trends in the be- 
haviour of the filling factors of emitting material with vari- 
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ous stellar parameters. The ability to reproduce these trends 
will then be a new constraint on theoretical studies of the 
structure and heating of late-type stars, including the mag- 
netic field geometry. 



8 CONCLUSIONS 

It has been shown that by assuming an energy balance be- 
tween radiation and thermal conduction in the upper tran- 
sition region, it is possible to use spatially unresolved stellar 
observations to estimate the fractional area filling factor of 
emitting material in stellar atmospheres. 

The technique has been applied to e Eri and it has been 
found that the area filling factor is ~ 0.2 in the mid tran- 
sition region and that it increases to ~ 1.0 in the corona. 
These filling factors are similar to those found from obser- 
vations of the solar transition region, and from models in 
which the transition region consists of material in magnetic 
funnels whose radii are greatest is in the upper transition re- 
gion/corona. Several potential sources of uncertainty in the 
calculations have been discussed and it is concluded that, 
although they are currently significant, they do not prevent 
a useful application of the methods adopted. Further obser- 
vations and calculations of atomic data should be able to 
reduce the current uncertainties. 
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